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ABSTRACT 

The linear stability of thin vertically-isothermal density-stratified Keplerian discs in 
toroidally-dominated magnetic fields is treated by asymptotic expansions in the small 
aspect ratio of the discs. The discs are found to be spectrally stable. The great vari- 
ety of possible initial conditions leads to three regimes of non-exponential growth of 
perturbations, which are classified according to different relative levels of the in-plane 
and axial perturbed velocities. The first two regimes of instability are characterized 
by the decoupling of the magneto-sonic (MS) and inertia- Coriolis (IC) modes, as well 
as by algebraic temporal growth of the perturbations, which are driven by either MS 
or IC modes (hereafter MS - and IC -regimes of instability, respectively). The third, 
mixed IC-MS regime of non-exponential, non-algebraic growth is due only to non- 
axisymmetric perturbations. The latter regime is characterized by high radial and 
azimuthal wavenumbers, and growth time of the order of tens of rotating periods. The 
mixed IC-MS regime most likely exhibits the maximal growth as compared with the 
IC- and MS- regimes. In the first two regimes of instability the compressible MS mode 
plays a principal role either as the driver of the growth or the driven growing mode, 
while the mixed IC-MS regime is described by the Bousinesq approximation for incom- 
pressible fluid. The latter is obtained as a natural limit of the expansion scheme. The 
presence of magnetic field in the mixed IC- MS regime may drastically increase the 
growth rates of the perturbations as compared with the pure hydrodynamic system. 
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1 INTRODUCTION. 

Instabilities in thin magnetized Keplerian discs are of principal interest in astrophysics. The numerous possible equilibrium 
magnetic configurations and initial levels of perturbations in real discs lead to a multiplicity of admissible regimes of insta- 
bilities, whose different aspects are shortly discussed b elow. Traditiona l ly, the results of the a nalytical study of the magneto- 
rotational instability (MRI) in infinitely long cylinders Velichoyl ll 19591); Chandrasechai 1 1960h ] are adopted in order to derive 
criteria for the normal-mode (spectral) stability in thin discs Balbus fc HawlevT i 1991 )]. These analytical results have been 
supported by num erical simulations, which d escribe the evolution of small perturbat io ns and their subsequent d evelopment 
to turbulence [e.g. iBrandenburg et al.1 |l995l h [rlawlev et all |l99rj ); iBodo et all l|200sh : iRegev fc Umurhanl |200sh ]. 

Adiabatic and vertically-isothermal disc mo dels. Two equations of s tate, adiabatic and vertically-isothermal ones, are 
commonly used for thin disc modeling [see e. g. iBalbus fc Hawlevl |l99lh ]. These two models correspond to quite different 
descriptions of the disc geometry: in the adiabatic model the disc has finite horizontal boundaries, while the horizontal 
boundaries of the vertically isothermal discs with the diffused mass density expand to infinity, and a finite effective disc 
thickness is properly defined. In the adiabatic disc model the hydrodynamic flow exhibits a relatively weak singularity of 
the basic functions at the horizontal edges of the unmagnetized discs [SMRRU2010], which becomes to be much stronger for 
magnetic field perturbations in toroidally dominated magnetic fields. The model of the vertically- isothermal disc is a fair 
approximation for thin discs under current consideration. It is characterized by the divergence of hydrodynamic velocities 



* E-mail: shtemler@bgu.ac.il; mond@bgu.ac.il; eliverts@bgu.ac.il 



2 Yuri M. Shtemler, Michael Mond, and Edward Liverts 



at infinity. Nonetheless, the perturbed mass flux, magnetic field and density vanish at infinity [see IShtemler et al.l (|2011h . 
hereafter SML2011 and the current study]. 

Incompressible- and compressible-fluid models. The assumption of incom pressibility of the pertur bed moti on is frequently 



emplo y ed in order to avo i d acou stic modes thus simplifying the system [e.g. iBalbus fc Hawlevl fl99lj), (19921. IColeman et al 
(|l995h . Curry fe Pudritz 1 199fj )]. Although in thin Keplerian discs the axial acoustic crossing time is of the order of the 



rotation period, in some of the cases incompressibility is a reasonable assumption. Indeed, for the poloidally dominated 
magnetic field configuration the acoustic waves decouple from the Alfven-Coriolis modes and, as was numerically established, 
converges to incompressible limit at least at sufficiently large plasma beta [see e.g. SML2011 and references therein]. For 
toroidally-dominated magnetic field, the current study demonstrates that the Boussinesq model for incompressible fluid is 
rigorously obtained in the limit of large plasma beta for the most dagenerous regime of instability [see Section 5 below]. 
However accounting for compressibility is crucial for other regimes of instability in toroidally-dominated magnetic field [see 
the current analysis in Sections 3 and 4]. 

The local plane-wave approximations and the asymptotic expansions in small disc 's aspect ratio. The simplest local approx- 
imation is the freezing equilibrium model that analyzes the shearing instability by fr eezing t he parameters of the equ i libriu m 
state in the Lagrangian cylindrical coordinates {r, 6, z} [e.g. Balbus fc Hawlev 1 199lh . ( 19921 . Ipapaloizou fc Teraueml l|l997h ]. 
Thin discs are thus modeled within the Boussinesq approximation for incompressible fluids ignoring the local structure ev- 
erywhere apart from the rotation radial shear effect that is the only inhomogeneity which is taken into account. Adopting 
Lagrangian shearing coordinates with O — 9 — fl(r)t (where f2(r) is the Keplerian rotation frequency), and applying the local 
approximation to the stability problem with respect to in general non-axisymmetric (spiral) perturbations allows solutions 
in the following form: /' ~ A(t) exp(zfc r r + ikO + ik z z). The effect of transition to the Lagrangian coordinates is reduced to 
replacement of a fixed radial wave number with a shearing time-dependent one, Ak r (t) = k r — ktdQ./dr. The time dependence 
of the radial wavenumber for non-axisymmetric modes does not allow a normal mode analysis, since it is inconsistent with 
Laplace transform in time. The resulting problem for the perturbations is solved by numerical integration in time of the 
corresponding ordinary differential equations with time-dependent coefficients. 

Summarizing the results of the freezing equilibrium modeling, axially-dominated magnetic field i s unstable on dynam i- 
cal (orbital) time of the disc rotation to both axisymmetric and non-axisymmetric perturbations [e.g. iTagger et al~1 \ 19921) ] . 
Whereas axisymmetric Alfven-Corioli s waves are expon e ntially unstable for a mixed toroidally-poloidal equilibrium magnetic 



configurations in high plasma beta Balbus fc Hawlev l|l99lh . Tagger et al. 1 1992h ] found that non-axisymmetric (spiral) 



perturbat ions are non- exponentia lly unstable for the poloidally-dominated magnetic field already at low plasma beta. Fur- 
thermore, Balbus fc Hawlev 1 19921 ) found that pure toroidal magnetic configurations are stable with respect to axisymmetric 
modes, but exhibit a non-exponential instability under non-axisymmetric perturbations, which occurs due to the time depen- 
dence of the radial wave number, Ak r (t). 

The shearing-box model is closely relate d to the freezing equ i libriu m approximation and is developed to study the 
shearing instability [see Umurhan et al. 1 20041 ) . Regev fc Umurhan 1 20081 ) and references therein]. The shearing-box model 
can be derived as a local expansion of the governing equations in the cylindrical reference frame and thus be expressed in 
Cartesian coordinates. It is widely accepted that the local shearing box approximation qualitatively reflects most necessary 
physics, such as differential rotation, compressibility etc. Transie nt growth of non-axisymmetric modes in pure toroidial 
magnetic field is exhibited within the shearing box approximation Brandenburg fc Dintrasl (200(3)]. Shearing-box modeling 
carried out for non-axisymmetric perturbations of a background pure toroidal m agnetic field within a rang e of resistivity 
and v is cosity, demonstrates inst ability of the linear Papaloizou fc Teraueml (1997)] and fully nonlinear system Fleming et al 



l200oj), Simon fc Hawlev ( 20091 ')]. Notwithstanding its success in predicting linear stability, the shearing box model has been 
recent ly criticized as to its viability in describing the ensuing nonlinear dynamical evolution of the system Umurhan et al. 

i 20oeh ]. 

Some results of the local modeling are still mainly qualitative due to ignoring the radial as well as axial structure of 
the disc. An approach alternative to the shearing box model introduced by iRegev (1983) is based on asymptotic expansions 
in a small aspect ratio, e, of thin discs. This approach allows to effectively describe the spatial inhomogeneity of the discs 
that makes it free from the locality constrains in both axial and radial variables. The thin disc approximation allows the 
adequate modeling of perturbations with the characteristic axial scale of the perturbations of the order of the disc thickness, 
in compliance with the commonly accepted view that turbulence is driven on the scales of the disc thickness. When true thin 
disc geometry is taken into account, both t he growth rates a nd number of spectrally unstab le MRI modes are greatly reduced 
and decreased with the disc aspect ratio Coppi fc Keve ] (|2003l ). lLiverts fc Mondl j2009h ]. These results for MRI strongly 
depend on the value of the plasma beta parameter [SML2011]. Thus, at large plasma betas starting from /? = iirP^/B^ ~ 5, 
the results of the boundary-value problem analysis tend to those obtained by plane-wave assumption in the axial direction. 
In that case, the discrete stability curve that describes the growth rate vs the e ffective axial wave numbe r for thin disc is well 
approximated by the corresponding continuous curve in the local modeling in IBalbus fc Hawlevl (|1991h . However, for beta 
values close to the minimal critical value the number of unstable modes is small and the disc stability properties significantly 
deviate from those predicted by the local model. Additionally, the thin disc approximation may be effectively combined with 
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the local approximations discussed above under additional restrictions on the plasma beta value (see Section 5 in the present 
analysis of discs in toroidally-dominated magnetic fields at large plasma beta). 

Non- exponential instabilities. The MRI is commonly thought as a proper source of the normal-mode (exponential) in- 
stability under either poloidally- or toroidally-dominated or mixed poloidal-toroidal magne tic fields. This assumptio n is sup- 
ported by analy tical r esults for infinitely long c ylinders or using th e relev ant local modeling Balbus fc Hawlevl Jl991 ), (1992), 
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(1996). 
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1 1997h . |Pessah fc Psaltisl |2005hJBrandenburg fe Dintras! l|2006h . lBegelman fc Pringlel l|2007l )]. For poloidally-dominated mag- 
netic fields the results of local modeling are consistent with those in thin disc approximation and exhibit the normal-mode 
instability (MRI) [SML2011], while thin discs in toroidally-dominated magnetic fields are rather exponentially stable. This ren- 
ders the system prone to non-exponential instability, which emerged as the dynamical response of thin discs either embedded 
in the toroidally-dominated magnetic field or even free from magnetic field at all. In general, two kinds of the non-exponential 
instability should be distinguished: algebraically and non-algebraically growth in time. It has been demonstrated that alge- 
braic instability is the primary source of pure hydrody nami c motion in thi n comp ressible adiabatic discs. For instance, this 

~~ i) and iRebusco et all (|2009h ] 



is establish ed both numer i cally 
in small e 



Umurhan et al. (200 



(2009)] and explicitly by asymptotic expansions 



Shtemler et all (|201ol ). hereafter SMRRU2010, where it is identified as the resonant as well as non-resonant in- 
teraction between inertia-Coriolis and sound waves]. Recently [SML2011] this approach has been generalized for Keplerian 
vertically-isothermal discs embedded in toroidally-dominated magnetic fields for one of the possible regimes (IC regime, see 
Section 4 below) of the instability investigated in the present study. It has been noted that such algebraic instability can 
lead to fruitful and phys i cally sound results in the context of the well-known scenario for sub-critical transition to turbulence 



Schmidt fc Henningsonl (2001)]. 



Poloidally - and toroidally-dominated magnetic fields. Although the magnetic field configuration in real discs is rather 
unknown, the poloidally dominated magnetic equilibrium configuration is commonly accepted in the initial stage of the disc 
rotation. MRI in an equilibrium under such circumstance s is generally thought as the instability that leads to sustainable 
turbulence in accretion discs [see iBalbus fc Hawlevl l|l99lh ]. However, according to observations and numerical simulations, 



the toroidal component of the mag netic field generated due t o the differential rota ti on may become comparable to or even 
dominate the poloid al one [see e.g. iBrandenburg et all j 19951 ). lHawlev et all |l996h . IStone et al. I jl996l ), lHawlev fc Krolik 



[ 2002h . |Proga"T[|2003h ]. This attracts attention to possible instabilities of Keplerian discs with a toroidally-dominated magnetic 
field. For axisymmetric perturbations, the existence of two regimes of instability distinguished by equilibrium magnetic 
configuration was demonstrated in thin disc approximation, namely either poloidally- or toroidally-dominated magnetic fields 
[SML2011]. In poloidally-dominated magnetic fields MRI (that is an axisymmetric normal-mode) was commonly found to 
be responsible for the most dangerous unstable mode, while under toroidally-dominated magnetic fields the non-exponential 
and non-axisymmetric governing modes of instability has been established in the current study. In SML2011 the toroidally- 
dominated magnetic configuration was treated for one of three admissible regimes considered in the current study: the 
non-exponential instability driven by either magnetosonic (MS) - or inertia-Coriolis (IC) -modes or mixed IC-MS regime. 

Main goals, assumptions and structure of the current paper. The current work is aimed to further elucidate the physical 
mechanisms responsible for the generation of the most dangerous modes in compressible density-stratified Keplerian discs 
in toroidal-dominated magnetic fields taking into account the true thin-disc geometry. Thus, the present study generalizes 
the stability analysis by SML2011 to non-axisymmetric perturbations and extends it to all possible regimes of instabilities. 
The analysis is carried out under the basic assumptions that a toroidally-dominated magnetic field is settled during the time 
preceding to the initial time in the stability analysis, and the disc instability operates on several orbital periods and on 
wavelengths of the order of or less than the effective disc thickness determined by the fluid compressibility (as the ratio of the 
sound velocity to the Keplerian frequency). 

The paper is organized as follows. The physical model for thin Keplerian discs that includes the dimensionless governing 
equations and their approximation to leading order in small aspect ratio for the steady-state disc are presented in the next 
Section. The governing equations and the results of analysis for perturbed thin discs in MS-regime of instability are presented 
in Section 3. Section 4 describes similar results for IC-regime of instability. Section 5 presents the stability analysis for the 
mixed IC- MS regime. The role of magnetic fields on the perturbation growth for three regimes is estimated in Section 6 by 
setting both equilibrium and perturbed magnetic fields to zero, and consideration of the pure hydrodynamic limit. Summary 
and discussion are presented in Section 7. 



2 THE PHYSICAL MODEL FOR THIN KEPLERIAN DISCS. 

The stability of radially and axially-stratified rotating plasmas in thin vertically-isothermal discs is considered. The discs 
are threaded by an equilibrium magnetic field whose toroidal component is larger than the other two components in a sense 
that will be defined below (the poloidally-dominated magnetic field is considered in details in SML2011). Viscosity, electrical 
resistivity, and radiation effects are ignored. 
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2.1 Governing equations. 

The dimensionless dynamical equations for vertically isothermal discs are the same as in SMRRU2010 and SML2011, and we 
briefly recall them here: 

DV_ 1 VP 1 j x B 

Dt ~ M| n /3M§ n ' [ ' 

^+V-(nV) = 0, (2) 

^+VxE = 0,V.B = 0, (3) 

E = V x B, (4) 

P = nT. (5) 

Here VP = c|Vn for vertically isothermal discs, and the dimensionless equilibrium sound speed is given by c| = dP/dn = 
T(r). Standard cylindrical coordinates {r, 0, z} are adopted throughout the paper; V is the plasma velocity; t is time; D/ Dt = 
d/dt + (V • V) is the material derivative; &(r, z) = — (r 2 + z 2 ) -1 / 2 is the gravitational potential due to the central body; B, 
j = V x B, and E are the magnetic field, current density and electric field, respectivelly; P = P e + Pi is the total plasma 
pressure; Pi = niTi axe the partial species pressures (I = e, i); T = T e — Ti is the plasma temperature; subscripts e and i 
denote electrons and ions, respectively. The positive direction of the z axis is chosen according to positive Keplerian rotation. 
The dimensionless coefficients Ms and j3 are the Mach number and the characteristic plasma beta, respectively: 

V* P* 

— ■ P = ^Tk 

cs* Bi 

The asterisk denotes the corresponding dimensional characteristic values defined along with the dimensionless variables in 
SML2011. In particular, V*, cs*, P* and B, are the characteristic velocity, sound speed, pressure and magnetic field, re- 
spectively. The characteristic radius value, r*, the disc semi-thickness, H(r*), the inverse Keplerian frequency, t* = f2 _1 (r*), 
and the background toroidal magnetic field, P* = Bg(r t ), are taken as the characteristic dimensional scales for radial and 
axial lengths, time and magnetic field, and from the basis for estimations of the variables for all admissible regimes under 
consideration (see Tables 1 and 2 below). Additionally, vanishing boundary conditions at infinity, z — ±oo, are adopted for 
the magnetic field and density perturbed about an equilibrium state. 

A common property of thin Keplerian discs is their highly compressible motion with large Mach numbers 
(2002)]. Furthermore, the disc aspect ratio e is a small parameter that is equal to the inverse Mach number: 



M S = -A p = 47r-|. (6) 



Frank et al 



The smallness of e means that dimensionless axial coordinate is also small, i.e. z/r* ~ e (\z\^H t ), and consequ ently the follow- 



ing re scaled quantities may be introduced in order to further apply the asymptotic expansions in e [similar to lShtemler et al 



1 20091) : SMRRU20fO; SML20U] 



C = --°, PTM = ^f4~e°, (8) 
where H(r) is the scaled semi-thickness of the disc; H(r*) = 1. 

2.2 Steady-state equilibrium configurations. 

It is first noted that the asymptotic expansion for the time-independent gravitational potential is given by: 
*(r,C) = Hr) + e 2 0(r,(), $(r) = ~, 4>{r, C) = ic 2 " 2 (r) + 0(e 2 ), r>i»e. (9) 
Substituting ((9]l into (0-© and setting to zero the partial derivatives with respect to time yield to leading order in e 
Vl = dMr± c s (r)dn_ ^(r,Q 

r dr ' n d( d( ' 1 ' 

Thus, the solution of l|10p is given by: 

Ve = e Ve(r) + O(e 2 ), n e°n = 6° N(r)u(jj), (11) 
and 

Ve(r)=rn(r), n(r) = r" 3/2 , u(r,) = exp(-r? 2 /2), r, = C/H(r), (12) 
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Table 1. S for the equilibrium gauge functions in H14I I. 

n V r Vg V z B r Bg B z 

> 2 > 2 > 2 1 



Table 2. S" for the perturbed gauge functions in IU4H and for the characteristic length scales of the pertur- 
bations in the radial, azimuthal and axial directions. 





n 


V r 


Vg 


Vz 


Bf 


Bg 


B z 


H T , Kg H z 


axisymmetric MS regimes 





2 


2 


1 


2 





3 


1 


non-axisymmetric MS regimes 





2 


2 


1 


2 





1 


1 


axisymmctric/non-axisymmetric IC regimes 











1 








1 


1 



non-axisymmetric mixed IC-MS regimes 0111000 1 1 



while other components of the steady-state velocity are of lower order in e: V r = o(e) and V z = o(e 2 ) [SML2011]. 

The equilibrium magnetic configuration considered below is characterized by scaling of the physical variables with e which 
generally may be written as f(r, = e s /(r, £). The magnetic field is assumed for simplicity to depend on the radius r only. 
It is also assumed that the equilibrium is characterized by a toroidally-dominated equilibrium magnetic field that is of the 
order e°, while the axial component of the magnetic field is of order e. These assumptions determine order in e of the rest of 
the physical variables [SML2011]: 

B r = o(e 2 ), Bg - e°Bg(r), B z - eB z (r),j r = o(e 2 ), jg - ejg(r) = -e^, j, <* e% = e° 1 (13) 

ar r ar 

All equilibrium variables are written in leading order in e, and depend on the radial variable only. The exceptions are 
the number density and the pressure that depend on the axial coordinate in a self-similar manner with radius-dependent 
amplitudes. In particular, the background toroidal and axial magnetic fields as well as the disc thickness and the amplitude 
factor, N(r), in the number density are arbitrary functions of the radial variable. Those functions specify the equilibrium state. 
Each equilibrium variable is characterized by a gauge function e s , where the parameters 5* for the case of toroidally-dominated 
magnetic field are presented in Table 1 according to (110[l - (|13[ l, 



2.3 Gauge functions for the perturbations. 

In the general unsteady case the dependent variables are scaled in e in the following way: 

f(rX,t) = e S f(r,C,) + e S 'f'(r,(,t). (14) 

Here / stands for any dependent variable, the bar and the prime denote equilibrium and perturbed variables, respectively. 
The perturbed part of each variable is characterized by a gauge function e s with some parameter 5*'. As shown below, the 
relevant regimes of the instability in such systems are the IC-, MS- and mixed IC- MS regimes, named so according to the 
mode driving the perturbation growth. The gauge functions for all regimes are summarized in Table 2 {j*c r , x-g and >c z are the 
characteristic length scales of the perturbations in the radial, azimuthal and axial directions). 



2.4 New dependent and independent variables. 

Before turning to the solution of the linearized stability system of equations it is noted that, guided by the steady-state 
solution and the special form of the dependence on the azimuthal variable, it is convenient to introduce the following new 
independent variables: 

r = ft(r)t, 6 = - ft(r)t, p = r -f- V = -J- (H(r) = ^M). (15) 

Jo H(r) H(r) Q(r) 

The derivatives in the new and old variables are related as follows: 

A - of )(— — - — — - 1 d JL - 1 [ d dlnti /d_ d \ ding d_i 

Here we keep without confusion the previous notations for Q and H, which now depend on the new radial variable p. It is 
also convenient to introduce the following new dependent variables scaled by radial variable dependent factors: 

V ri B' 

v(p,7J,t) = — rr, v(j>,t),t) = - , b(p,r?,r) = . (16) 

cs{r) N(r)v(ri) Bg(r) 
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In addition, the following functions of the new radial variable p are defined: 

2m_ a N{p)c%{p) . _ B e (p) _ r (f r (/A _ <jg(gK 

B z (p) H(p) Q(p) 

= , . dlnif - , . dlnfi - , , dlniV - , . dlnB| = , . d\nB z . , 

D H (p) = — Dn{p) = Dsip) = — D Be ( P ) = D bM = (17) 

Since the disc thickness is widely undetermined, and the above equations are significantly simplified in the case of disc of 
constant thickness, this will be adopted everywhere below 

H{p) = l, D B (p) = 0, pD Q {p) = -3/2, f(p) = r = p. (18) 



3 THE MS-REGIME OF INSTABILITY. 

3.1 Governing equations for the MS-regime. 

We start by substituting the decomposition of the total disturbed variables (|14p into the MHD system ([1} - (0 , and linearizing 
the resulting equations about the steady-state equilibrium solution (|10p - (|13p . By using the estimations for the equilibrium 
variables and the perturbed variables for the MS-regime (see Tables 1 and 2, respectively), the linearized momentum and 
mass balance equations are given to leading order in e in the following form: 

^ + nw ^^ w ._ 4w - ( . ) + _L_ s , (r) M, 

Maxwell's equations are now conveniently analyzed separately for the axisymmetric and non-axisymmetric regimes. 
Axisymmetric MS-regime. 

gt^'y-,, (25) 

at r or 

Note that the equations for both components of the perturbed poloidal magnetic field are decoupled from the rest of the 
system of equations, and may be dropped out from further consideration. In that case Eqs. (|23[) and (|25|) may be reduced to 
a single equation for the magnetic flux function, which guarantees the divergence-free condition VB' = 0: 

X = -±%; B^- 9 -f. (26) 
r oQ r or 



Although this means that Eq. ()23[) for B' r can be dropped, and Eq. (I25[) for B' z is sufficient for the further analysis, it will be 
convenient to calculate B' r and B' z directly from Eqs. ()23p and (12511 . 

Non-axisymmetric MS-regime. 

-8r +n ' r >-85"- B, < r )7-s5-- a(r) "8r = ■ ( ' 

These equations for thin discs embedded in toroidally-dominated magnetic fields (B r = 0, Bg ~ e°, Bz ~ e) are quite different 
from their counterparts in the axisymmetric MC-regime. This is so not only due to the presence of terms with d/d9, but also 
due to the appearance of B' z in the equation for B'g. Indeed, the e-scaling of B' z is different in those two cases: eB' z ~ e in the 
non-axisymmetric case, while tB' z ~ e 3 in the axisymmetric case (see Table 2). 
Applying the operator (V-) to the left-side of ([27|) - (f29|) yields 



Non- exponential growth in thin Keplerian discs under toroidally- dominated magnetic fields 7 



- fi(r)^^ + fl(r)V ■ B' = 0{e\ (30) 

where 

Thus, the divergence-free condition, V • B' = 0, satisfied at the initial instant, is guaranteed to be satisfied for all subsequent 
times in consequence of the Maxwell equations (|27|) - (|29[) . Consequently, finding from the divergence- free condition dB' z /dC, = 
—r~ 1 dB' g /d8, and substituting it into equation ()28[) the latter may be rewritten in the standard form with the first two terms 
as material derivative of B' g : 

It is first noted that the radial component of the perturbed magnetic field e 2 B' r ~ e 2 is small compared with other 
components of the perturbed magnetic field: e°B' g ~ e° and eB' z ~ e. Moreover, the input of B' z into the divergence-free 
condition V ■ B' = is of leading order in e, while the corresponding input of B' r is negligibly small: 

2 1 d(rB' r ) ldB' e dB' z 

Furthermore, B' r can be set equal to zero without loss of generality 

B' r = 0. (34) 

Consequently Eq. (|34|) substitutes (|27p for non-axysimmetric perturbations, and the magnetic field is two-dimensional to 
leading order in e, B' = {0, e 2 B' g , eB' z }. As a result, the magnetic flux function ^f' (different from that introduced in (I26|) for 
the axisymmetric case) can be used so that the divergence-free condition is identically satisfied: 

B '-7W Be -~^t> (35) 

and (|29[) for B' z can be dropped from further consideration. Again, however, it is convenient to find B' z and B'g directly from 
([29]) and (32). 

Examining the resulting system of equations, namely, (|19[) - (|22[ ) and either (I23I) -(I25 [ ) (in the axisymmetric case) or ((27])- 
(1291) (in the non-axisymmetric case), it is seen that the entire system decouples into three subsystems: (i) for the MS mode 
that consists of equations for the perturbed density, n', axial velocity, V z , and toroidal magnetic field, B' e ; (ii) for the IC 
mode that consists of equations for the perturbed in-plane velocities l/ r ' and B' ; and (iii) for the perturbed poloidal magnetic 
fields, B' r and B' z . The systems of equations for the MS and IC modes have an invariant form for both axisymmetric and 
non-axisymmetric perturbations, as distinct from equations for perturbed poloidal magnetic field which should be described 
separately for axisymmetric and non-axisymmetric modes. Thus, the above mentioned subsystem (iii) for the perturbed 
poloidal magnetic fields can be dropped in the non-axisymmetric case, since B' r = and B' z can be expressed through the 
magnetic flux function (|35[) or, alternatively, through B'g, while in the axisymmetric case, introducing the magnetic flux 
function (|26|) allows to reduce the subsystem to only one component of the perturbed poloidal magnetic field, e.g. B' z . Let us 
now rewrite the governing system for the MS mode in terms of the new variables (I15|) - (|18[ ), 

The MS mode of the MS-regime. The equations for the in-plane components of the magnetic field decouple from the above 
system of governing equations, while the perturbed toroidal magnetic field, B' B , together with the equations for the perturbed 
density and axial velocity, n! and V z , respectively, decouple from the equations for the IC mode, i.e. for V r ' and Vg. 

Thus, in the new variables (|15p ~ (|18p the MS subsystem of equations yields for the perturbed axial velocity, density and 
toroidal magnetic field, u, v z and bg: 

8u + 1 dju&M 
or 1/(77) or\ 

^T + ^=0. (38) 
or drj 

Note that all equations for MS mode have an invariant form for both axisymmetric and non-axisymmetric modes, and are 
also free from radial derivatives. 

The IC mode of the MS-regime. Similarly, the following subsystem for perturbed in-plane velocities describes the IC 
modes: 

dv r „ =, , . ( dv 1 1 dbg 1 d\v(rj)vz\ 1 11 dv z 1 1 , „ T .„„. 

or I dO j3g{p) u{rj) d& 1/(77) dr] I Pe{p) v(r)) dr) Pe{p)v(rf) 
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dve 1 = 1 dv , _ 1_ 1 dfo , 4f) x 

5r + 2 Ur P dQ + S{ P )pe{p) Hv) d V ' 1 ' 

where 

<%e rl - . , 2i, ar dv 1 1 1 rl - . . 2i 

Lobe s i^™ w+ pJ 69 - ^ s -^+2 m m ^ m (p)+ pJ " 

with v, v z and &e in the right-hand sides of (l39|) - (f40|) determined from (|36|) - (f38)l . As it is evident from ({39]) - (f40|). the MS 
mode acts as a driving force (hence the term MS regime. Note that the first term in the right-hand side of (|39|l is proportional 
to time. It arises from the radial derivatives of the perturbed density and toroidal magnetic field due to the radial shear of 
the Keplerian angular velocity (i.e. Dci(p) ^ 0). 

The poloidal magnetic field in the MS-regime. As mentioned above the equations for the perturbed poloidal magnetic field 
are described separately for the axisymmetric and non-axisymmetric modes by Eqs. (I23[) . (|25[1 and (J29J) , (|34p . respectively: 

(i) The axisymmetric case. 

db r 1 dv r 



dr S(p) dn 



(41) 



db z 2 Dl(p) 1 d[v(Ti)v M ] Dn(p) 1 r 1 0», , . - , ,, .. , ■ ,. 

T - — r^r \ -= -, : -j. h lv(r\)ve + - Lobe + iVoH — - Mov r , (42) 



9r S(p) dr, S( P ) Hv)^o(p) dv Mp) J SQ>) 

where 

A/ « r = ^+\Dn(p)+~DBz(pj\ v r+-Vr- 

op L 2 J p 

(ij) TTie non-axisymmetric case. 

b r = 0, (43) 
<% z _ 1 dv z 

~d7 ~ p ~Be' 

Thus, for the axisymmetric modes the amplitudes of the perturbed poloidal magnetic field, b z , grows quadratically with time, 
while in the non-axisymmetric case it is constant. 



(44) 



3.2 Normal-mode stability for the MS-regime. 

We start by representing the perturbations as follows: 

f(p, ©, V, t) = expf-iAr + imO]f(p, rf) + c.c. (45) 

Here / stands for any perturbed variable, A is the eigenvalue, and m = 0, ±1, ±2, ... is the azimuthal wavenumber. Since the 
resulting equations which describe the temporal growth of the perturbations do not contain radial derivatives, the normal- mode 
form (|45l) for the characteristic length scale x r ~ e° is used without assuming an exponential form in the radial direction. 

The MS mode of the MS-regime. Substituting the equilibrium number density v(r)) = exp(—r) 2 /2) and the ansatz (|45[1 
into the subsystem ({36} - (J38j) for the MS mode, the following single ordinar differential equation is obtained for the toroidal 
magnetic field bg: 

d 2 bg Pevjrj) - I dbg ^ PeHv) 



, V 7^'_{ | - ~£ + (A - 2) - T-r b S = 0, b e =0, at r, = ±oo, (46) 



while the axial velocity and density, v z and v, can be easily expressed through bg: 

i) z = i\ I bgdr/, v — bg — tj / bgdn. (47) 



These equations for the MS mode in the MS-regime are similar to those in SML2011 for toroidally dominated magnetic fields. 
Summarizing those results we represent the approximate solution of (|46[l obtained by applying the WKB approximation: 



bo = Q( V )eM[ f*d V ), p= ^f_^ 1 (48) 

The calculations are simplified if it is additionally assumed that fig 2> 1. This yields 

W f W f 

exp[±i / x(rf)drj\ and Q = — nxp[— / | x(r\) \ drj\. (49) 



Approximate expressions for p,(rf) and >c(r\) to leading order in (3g are 
H(V)M~V, * 2 (v)*-~(ri 2 -rit), (v. « 2 | A |). 
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Equations (|49[) describe the oscillating solutions for *c 2 (r]) > 0, | 77 |< 77* and monotonically decreasing for x 2 {ri) < 0, 
77* <[ 77 |< 00, respectively, where the turning points 77 — ±77, are the zeros of >t{rf). 

To completely determine the solution, the matching conditions betweenthe inner (| 77 |< 77*) and outer (| 77 |> 77,) regions 
are applie d in order to obtain the co efficients in (|48|) and the eigenvalue equation. The result is the following Bohr-Zommerfeld 
Landau fc Lifshitsl (|l997T )] that determines the dispersion relation for the eigenvalues A 



condition 



+1* 1 r - 

x(r\)dr) = n(k + -), and A = X MS = ±Vk + 1 w ±Vfe, (50) 
n* 2 

where k is the number of zeros of the solution for be, and (|50[) is valid for sufficiently large values of A and k, as required in 
the WKB approximation. 

The IC mode of the MS-regime. For the non-resonant case with non-equal frequencies of IC and MS modes, the perturbed 
axial velocity, torodial component of the magnetic field and density are zero, i.e. 

v z = v = b e = 0. (51) 

Therefore, the perturbed in-plane velocity is described by a set of equations that governs the dynamics of the IC modes. 
Substituting the ansatz flM) and (O into (O - @3 yields [SML2011]: 

— i\v r — 2ve = 0, ~i\ve + ^v r = 0. (52) 
The corresponding dispersion relation is therefore given by 

A = \ic — ±1, (53) 

which represents stable epicyclical oscillations in the disc plane. As the axial and radial coordinates play the role of passive 
parameters, the eigenfunctions of the IC modes are determined up to arbitrary amplitude A(p, 77) from (|52l) - (|53l ): 

ve=T^iVr = =F~iA(p,ri), A(p, v ) = F{p)G{ v ). (54) 

Since any special form of the function A(p, rj) represents some given set of initial conditions, the self-similar separable form 
of the planar velocities with the arbitrary functions F(p) and G{rf) is considered as an example. 



3.3 Algebraic temporal growth of perturbations in the MS -regime. 

A detailed account of resonant as well as non-resonant excitation that lead to quadratic and linear temporal growth, respec- 
tively, is given for pure hydrodynamic systems [SMRRU2010 and SML2011]. Here, however, the eigenvalues Xic = ±1 01 the IC 
normal modes in ([S3J do not coincide with those for the MS modes in (f5H)l . A A /s ~ ±vfc, in the large k approximation adopted, 
Xms 7^ Xic, the resonance condition is not satisfied. Focusing therefore on the IC modes that are non-resonantly driven by 
MS modes {be,v z ,C / }, and consequently on the linear temporal growth of the amplitude, the variables that characterize the 
IC-modes are given by: 

{v r , ve} = r{'fi* 1) ,75^ 1) }(p, 77) exp(-iA A / S r + imQ) + c.c, {b r , b z } = {r 3 bi j) , r l b z l) }{p, 77) exp(-iA A ,/sT + im&) + c.c, (55) 

where m = 0,±1,±2, ... is the azimuthal wavenumber; Vr , Vg 1 ^ and W\ VP are amplitudes that depend on p and 77, while 
the parameters j and I will be specified below; the radial magnetic field is either zero for non-axisymmetric case or, as 
for axisymmetric case, can be found from Eq. (|41|l . Substituting (|55|l into the system (|39p - (|40p for v r and vq, and keeping 
only terms of the highest powers in r (the rest of the terms in the solutions are dropped here for brevity), yield for both 
axisymmetric and non-axisymmetric cases: 

~~[xms = 2v ° M = >^i Daip) v m [ v+ r 9 ^n) be \ + W)^r$' (56) 

where v, v z and be are given by (47) and (48); the terms proportional to azimuthal wavenumber m and shear coefficient Dn(p) 
arise due to the joint non-axisymmetric and shear effects, while the rest of the terms proportional only to Dq (p) are entirely 
due to shear. Substituting the normal-mode ansatz (|45|) into (|41[) -(|42 | ) for poloidal magnetic fields b r , b z in the axisymmetric 
case (m = 0) and keeping only terms of the highest power in r yield 

eW=<^^, 8«m = ^-^ ^ »0M. (57) 

S(p) or) iXms S(p) u{rj) drj 

Similarly, ([43j) and ()44[) for b r , b z in the non-axisymmetric case (m 7^ 0) yield 

6( 0) (p,77) = 0, bf\p,rf) = (58) 

Aa/s p 
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Thus, for both axisymmetric and non-axisymmetric cases the equilbrium toroidal magnetic field gives rise to algebraic insta- 
bility that is driven by the MS modes. In the axisymmetric case the in-plane perturbed velocities and the perturbed radial 
magnetic field grow linearly in time, while the perturbed axial magnetic field grows quadraticaly in time. In turn, in the 
non-axisymmetric case there is a linear temporal growth of the in-plane perturbed velocities without algebraic growth in the 
perturbed magnetic field. The interesting fact is that the non-axisymmetric modes of the MS-regime of instability do not go 
over to the axisymmetric value at small azimuthal wavenumbers. This occurs due to nonuniform behavior of the perturbed 
axial magnetic field in e that is ~ e and ~ e in the non-axisymmetric and axisymmetric cases, respectively. 



4 THE IC-REGIME OF INSTABILITY. 

4.1 Governing equations for the IC -regime. 

We again start by linearizing the MHD Eqs. JT]) - <[6j) about the steady-state equilibrium solution (110p - (|13[ ) and substituting 
the appropriate expansion of the total disturbed variables (|14|l for the equilibrium and perturbed variables for the IC-regime 
(see Tables 1 and 2, respectively). Keeping terms of leading order in e the system of equations is decoupled into the following 
subsystems reformulated, as for the MS-regime, in terms of the variables defined in (|15[) - (|18l) : 
The IC mode of the IC-regime. 

^--2v e = 0, (59) 

The MS mode of the IC-regime. 

dvz dv 1 1 db e = . . 

Or + dr, + Mp) Hv) d V 1 ' 

dv 1 d\v(rf)V;,] = , . \dv r n 1 \dvg i dv r ,- , . - , .. , 

d~r + W)^r = TDn{p) la£ - H ~ - P lid + Vr ] ~ a? " [Dn(p) + DN(p)]Vr ' (62) 

dbg dv z = . .\dv r i 3 1 dve dv r r = , . 1= , , . 

The poloidal magnetic field in the IC-regime. 

8br__'\_dVr_ 1 dv r 

dr ~ p 96 S(p) d V ' ( ' 
db z D n {p) \dv r 1 1 dv z 1 [1 dv g dv r / - 1 - , 1 

-d^ = r -JipJ Ido 2ve \ + ~p ae - W) U> do + ~p - \ Da{p) + 2 Dbz{p) + p) Vr \ ■ (65) 

Again, if the divergence-free condition, V • B' = 0, is satisfied at the initial instant, it will be fulfilled for all subsequent times. 
Additionally, Eq. (I63|l for the perturbed toroidal magnetic field was derived by employing the divergence-free condition of the 
magnetic field. The system of Eqs. (159[) - (|65[) is subject to the vanishing boundary conditions at infinity for the perturbed 
density and magnetic field. 



4.2 Normal-mode stability for the IC-regime. 

In order to describe the IC -regime, we start again by representing the perturbations in the normal-mode form for any 
perturbed variable /, with the complex eigenvalue A and the azimuthal wavenumber m — 0, ±1, ±2, ... 

f(p, 9, V, T ) = exp[-iAr + imQ]f(p, n) + c.c. (66) 

Substituting that ansatz into the governing equations reveals that the spectral properties of the IC-regime are identical to 
those that characterize the MS-regime. Indeed, the subsystems (|59[) - (|60|) and (|61|l - (|63|l for the IC and MS modes of the 
IC-regime are described by (|52p - (|54|l and (146|l -(|50 [ l for MS-regime in Section 3. However, as will be seen in the next sub- 
section, while in the IC regime the MS modes drive the IC waves to non-exponential growth, the opposite happens in the 
MS regime. As in Section 3, the resulting equations which describe the temporal growth of the perturbations don't contain 
radial derivatives, the normal-mode form (66) is used without assuming an exponential form in the radial direction with the 
corresponding characteristic scale *c r ~ e°. 



N on- exponential growth in thin Keplerian discs under toroidally- dominated magnetic fields 11 



4.3 Algebraic temporal growth of perturbations in the IC-regime. 

Assuming that the MS modes are non-resonantly driven by IC- modes (that is determined by (|52[) - (|54[) ~) and substituting 
(1671) into in the right-hand sides of (I61|) - (|63[) for the MS mode and (|64l> - H65[) for the poloidal magnetic fields {b r ,b z }, yields 
(keeping for brevity only linear in r terms in the right-hand sides (|59[) - (|65[) ): 



{v z ,v,b r ,bg,b z } = t{v z 1 \ u {1 \b r 1 \b e \ 1 \b z 1) )}(p, ri) exp(-iX Ic r + imO) + c.c. (67) 

Assuming again non-resonant conditions and using Xic — ±1, yields after cumbersome calculations the following partial 
solution that is driven by the IC mode in the right-hand sides (159[l - (|65[) : 



6}. j = 0, by = - — - 7 ' y / (imt;r-2u 9 )di/(7?), ^ = -- — g^-r (imv r -2v e ), 

v z 1} = 5n(p) / (tmfir ~ 2fi e )dr? + *A /C / S< 1 >*7, » W = ^ + ^^T^T^- (68) 
J J iA/c f(r?) dry 

Expressions (|68[) determine the perturbed magnetic field, density, and axial velocity, through the radial and toroidal velocities, 
(imv r , vg), which are the eigenfunctions for the IC mode introduced in (|54f) . arbitrary functions, F(p) and G(p), should be 
additionally specified in order to satisfy the corresponding boundary conditions. 



5 MIXED IC-MS REGIME WITH HIGH AZIMUTH AL AND RADIAL WAVENUMBERS. 
5.1 General magnetohydrodynamic system. 

Within each of the two regimes discussed above, namely the IC and MS regimes, the two modes of wave propagation (IC 
and MS waves) are practically decoupled, as they may exist independently of each other. In addition, within each of the 
regimes one of the modes may also drive the other. To avoid such mode decoupling, the effect of large values k r , kg ~ k z ~ e 
is examined in accordance with the principle of least possible degeneracy of the problem (where k = {k r ,ke,k z } is the wave 
vector, k r — 2n/ ' n r , kg — 2n / Kg, k z — 2-k/k z , see the last row of Table 2). As is subsequently shown for perturbations of such 
scale lengths, both modes are coupled, and the disc exhibits a strong instability which is absent within IC and MS regimes 
discussed above. As the radial as well as the azimuthal scale lengths of the perturbations are of order e, the perturbations are 
represented in the following normal form in the axial and the radial variables: 

f{p, 6, r], t) = exp[i(k p p + me6 + k r} rf)]f{p, r), t) + c.c, (69) 

where to leading order in e, f(p, rj, r) parametrically depends on p, and is independent of r\ (i.e. the disc structure in r\ may 
be ignored), 

k p = ek r ~ e°, fee — efce = me//) ~ e°, (me — era ~ e°), k v = ek z ~ e°. 
Employing now the expansions described in the last row of Table 2, yields: 

^ - 2vg = -iAk p (r)i> + -L -L{[Ak 2 p (r) + k%]b r + Ak p (T)k v b z }, (70) 
or fc e p e 

+ = -ikev, (71) 

^ = ^ik z i> + -^--L[(k 2 e + kl)b z + Ak p (r)k v b r ], (72) 
or fee fig 

Qv 

— = — ikzVz — iAk p (T)v r — ik&v e , (73) 

OT 

^ = ik e v r , (74) 
or 

— — ik z v z — iAk p (j)v r — ~8 r , (75) 
^ = ifceu z , (76) 

OT 

where Afc p (r) = k p — TpDn{p)ke{p). Equations (J72j| - (JTHJl imply that 
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Hence, the divergence-free condition is satisfied at all times by choosing a proper initial condition for the magnetic field. 

Assuming now large plasma beta, and applying the asymptotic principle of least possible degeneracy of the problem in 
the limit fig S> 1 leads to the following estimations: 



fee ~ 


• Pi are 


introduced 
















kp — 


kp 
fee ' 


I k& - 1 
fee = -r- = 1, 
fee 


fe?7 — 


k-q 
fee ' 


Afe p (T) 


Afc p (r) 
fee 


■' Afc(T) " fee 




= & 
fee 


Vp = 


vk%, 


b r — b r , be 


= 6e, 


k = 


bg, V r 


= kev r , 


«e = keve, v z = 


kev z , 





fc p ~ fc© ~ fc„ ~ Afc p (r) ~ yjp g , V^PgX ~/3g%. (78) 
Choosing now for further convenience fee as the basic parameter instead of fie , the new scaled variables of order of unity in 

(79) 

(80) 
where 

Afc 2 (r) = Afc^(r) + l + fc^, Afc p (r) = fc p -pL>n(p)r = fc p + ^r. (81) 
This yields the following equations with coefficients of order of unity in the plasma beta parameter: 

- 2ivg = Afc p (r)*> - -L{[Ak 2 p {r) + l]b r + Ak p (T)k v b z }, (82) 

OT fig 

dive 1.- - ,™> 

-J. h TTlVr = v, (83) 

OT Z 

^ = fc^j? _ -L[(k 2 v + 1% + Ak p {r)k v b r ], (84) 
or pe 

iAkp(T)v r + ivg + ik v v z = 0, (85) 

dhr ;~ fans 
— — = iv r , (86) 

OT 

^ = -ik v v z - iAkp(r)vr - |&r, (87) 



db z 

TJ = IVz 
OT 



(88) 



Here the parametric dependence on p is omitted for brevity; the time derivative of the perturbed density, D, in (|85(l is neglected 
due to its smallness at large plasma beta fig 3> 1, while the steady-state (|85ft is reduced to the divergence- free condition for 
the plasma velocity. It is interesting to note that condition (|85|) . which is tantamount to the Boussinesq approximation, is not 
assumed, but rigorously obtained from the expansion scheme in small e and 1/fig ^> 1. It can be directly verified using (18611 
- (1881) that the divergence- free condition (|85[) for the velocity v is satisfied identically due to the divergence- free condition for 
the magnetic field, rewritten as 

iAfc p (r)6 r + ibg + ik v b z = 0. (89) 

Equations (|82|) -(JHHJ) may be now reduced to the following coupled system of equations for the poloidal magnetic field 
and the density perturbations: 

ft = -2k v ^ - 2Afc p (r)^ + Ak p {T)D - ^{[AP p (t) + l)b r + Ak p {r)k n b z }, (90) 

OT CJT C/T /3q 

^-f = ~k v u - -^[(1 + k\)b z + Ak p (r)k v br], (91) 

OT p g 

1 [2Afcg(r) - A + 2 fc ^ A W g» + tAk p (T)b r + tk n b z . (92) 



Afc 2 (r) P 9r Afc 2 (r) <9r /3 e /3, 

Substituting v from (|92[l into (|90|l - (|91|l finally yields a system of two coupled second order differential equations: 

ft + (3 + 2^)4M| +2^1+^ + ^ = o, (93) 
<9r 2 ^Afc 2 (r) <9r ^Afc 2 (r) 8t p g ' X ' 

ft - fe f^)- 1 ^ -2%mi ^ + U. = 0. (94) 
dT 2 71 Afc 2 (r) 9t n Ak*(T)dr fa 

After b r and b z are found from (|93[) - (|94|) . the density is determined by (|92|) . and the rest of the variables v r , v z and vg, bg 



can be determined from the (|851) - (|86|) an d (L5BJ) ~ (|89l) . respectively 
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Figure 1. Amplitudes of perturbed magnetic field, velocity and density vs orbital time in the mixed IC- MS regime of instability 
for typical parameters k p = 0.5, k v = 1.5, j3g = 0.5 and partial initial data b r (0) = 0, iv r (0) = db r /dr(0) = 0, b z (0) = 1, iv z (0) = 
db z /dr(0) = —1. (a) poloidal magnetic fields, b r (solid line) and b z (dashed line); (b) toroidal magnetic field, bg; (c) poloidal velocities, 
iv r (solid line) and iv z (dashed line); (d) toroidal velocity, vg; (e) density, v. 
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mm 

(6> (c) 

Figure 2. The same as in Fig. 1 for k p = —0.5, = 1.5, 0$ = 0.5 






IVg 



-Afe p (r) 



db r 



k dbz 



ibg = — Afc p (r)b r 



(95) 



db r <9b z 
iv r = — , m; z = — — 
or or 

Numerical solution of (|92[) - (|95p for amplitudes of the perturbed magnetic field, velocity and density vs the orbital time 
t in the mixed fC- MS (oscillatory) regime of instability are presented in Figs. 1-4 for typical values of the scaled (with 
the azimuthal wave number) parameters, k p , k v , fig. It is noted that without loss of generality k v may be assumed to be 
positive since changing its sign is equivalent to sign changing of either b r or b z . Figures 1-4 demonstrate that the axial 
magnetic component b z dominates the radial one b r for all values of the parameters. Moreover, a significant amplitude growth 
of the axial magnetic field b z may be achieved after tens of orbita l times. It should be not ed that the results obtained for 
the mixed IC-MS regime of instability are quite similar to those in iBalbus fc Hawlev (1992). Such a growth of b z may have 
a non-monotonic character which is evident by its decrease with further increasing of time, i.e. the instability may be of 
transient nature. Figures 1-4 also exhibit strong dependence of the magnetic field, velocity and density on the values of the 
parameters. It is noted that increasing fig and k v has a dramatic effect on amplification of the initial perturbations. Thus, 
the pure toroidal magnetic field configuration is susceptible to high beta, high small axial wave length, and low azimuthal 



wave numbers instabilities. In particular, a significant amplitude growth of the perturbations by factors ~ 10 4 
achieved after tens of orbital times (see Figures 4 and 5). 



10 may be 



6 THE PURE HYDRODYNAMIC LIMIT. 

To elucidate the role of the magnetic fields in the perturbations growth, consider the pure hydrodynamic limit by setting both 
equilibrium and perturbed magnetic fields to zero, while maintaining the same assumptions and notations in all above three 
regimes. In that case the regimes of instability may be classified by a quite similar way as for the general MHD case. There 
are three following regimes: the MS regime (degenerated in the pure hydrodynamic limit to the acoustic regime) and the IC 
regime, which are characterized by algebraic temporal growth of the perturbations, as well as the mixed IC-acoustic regime 
for which the non-algebraic instability has been established for non-axisymmetric perturbations. Below, we restrict ourselves 
the mixed IC-acoustic regime as the most degenerous one, while the first two regimes are not described here for brevity. 
The mixed IC-acoustic regime. In a similar way to the mixed IC-MS-regime the pure hydrodynamic limit is as follows: 




(a) W 6A CS) <e) 



Figure 3. The same as in Fig. 1 for k p = 3.5, k v = 12, /3g = 3.5 . 
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Figure 4. The same as in Fig. 1 for k p = —3.5, k v = 12, 
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Figure 5. Amplitudes of perturbed velocity and density vs orbital time in the mixed IC- MS pure hydrodynamic regime of instability 



for typical parameters k p 



0.5, k v 



1.5, 



0.5 and partial initial data iv r (0) = 0, ii> z (0) = —1 (or, equivalently, i>e(0) = k^). (a) 



poloidal velocities, iv r (solid line) and iv z (dashed line); (b) toroidal velocity, vg; (c) density, v. 



dr 2 Vr Afc 2 (r) 



(v r + 2vg), 



Ak 2 (r 



■(v r + 2ve). 



(97) 
(98) 

Equation (|98|l expresses the divergence-free condition for the perturbed velocity. The results of the numerical solutions of Eqs. 
H96 [l -H98 [ ) are presented in Figs. 5-8 for the same values of the parameters as for the MHD system in Figs. 1-4, respectively. 

Comparing equations for the MS- and IC-regimes in the general MHD systems with their counterparts in the pure 
hydrodynamic limits, it can be shown that at large plasma beta the magnetic field has no significant influence on the 
linear temporal growth of the in-plane velocities. In contrast, the behavior of the pure hydrodynamic and MHD systems 
is quite different for the mixed IC- MS regime, where the presence of magnetic fields results in a dramatic increase of the 
perturbation growth rates. It s hould be noted that the pure hydrodynamic li mit presented in this section generalizes previous 
hydrodynamical results such as lUmurhan et al. | j200rj ), lRebusco et all (|2009r i and SML2011. The latter focused only on those 
perturbations which belong under the current classification to the IC regime (i.e. are scaled according to the fourth row 
in Table 2). Indee d, typical amplificatio n of such perturbations do es not exceed 1-2 orders of magnitude after about 100 
rotating times [see lUmurhan et al. (2006) and lRebusco et al.l (|2009l )]. This is in contrast to the much higher amplification of 
perturbations that belong to the mixed regime under the influence of a toroidal magnetic field, which may exceed six orders 
of magnitude (compare Figs. 4 and 8). 



7 SUMMARY AND DISCUSSION. 

The present study concerns the stability of thin density-stratified vertically-isothermal Keplerian discs in a toroidally- 
dominated magnetic field. The problem is treated by asymptotic expansions in small aspect ratio of the disc. The vertically- 
isothermal thin disc approximation is found as the most suitable model to stability study of Keplerian discs due to its relative 
simplicity and adequacy to true disc geometry. The model may be effectively combined with the local approximations (as it is 
done here for the mixed IC-MS regime). The toroidally-dominated magnetic configurations are found to be stable with respect 
to normal-mode perturbations. Instead, three non-exponential regimes of temporal growth of perturbations governed by the 
rotation shear have been identified depending on the initial level of perturbed velocities. First two regimes (magneto-sonic 
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Figure 6. The same as in Fig. 5 for k p = —0.5, k v = 1.5, /3g = 0.5 
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Figure 7. The same as in Fig. 5 for k p = 3.5, k v = 12, /9e = 3.5 
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Figure 8. The same as in Fig. 5 for k p = —3.5, k v = 12, fig = 3.5 



(MS) and inertia-coriolis (IC) regimes are driven by MS- and IC-modes, and decoupled into in-plane and normal-plane modes 
with small in-plane components of velocity compared with axial velocity component and contrariwise, respectively. The MS- 
and IC- regimes support the perturbations with algebraically growing amplitudes for both axisymmetric and non-axisymmetric 
modes. The third, mixed IC-MS regime is pure non-axisymmetric, operates on high radial and azimuthal wavenumbers, and 
corresponds to the least possible degeneration of the problem with comparable in-plane and axial perturbed velocities. This 
regime most likely corresponds to maximal growing perturbations, which demonstrate significant growth rates with a growth 
time measure d in tens of orbital period s. The corresponding time-dependences of the magnetic fields are qualitatively sim- 
ilar to that in iBalbus fc Hawlevl 1 19921 ). Moreover, although the compressibility effects are taken into account consequently 



through th e current analysis, the div ergence- free condition for the Bousinesq approximation is valid in both models: as a'priori 
accepted in lBalbus fc Hawlev (1992) modeling, or as the corollary of the high plasma beta approximation in the present anal- 
ysis. In the first two regimes of instability the compressible MS mode plays a principal role either as the driver of the growth 
or the driven growing mode, while the mixed IC-MS regime is described by the Bousinesq approximation for incompressible 
fluid that is obtained as a natural limit of the expansion scheme. Additionally, as distinct from IC and MS regimes, behavior 
of the pure hydrodynamic and MHD systems for the mixed IC- MS regime, the presence of magnetic fields may drastically 
increase growth rates of perturbations. 
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